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The early observation at HERA of an excess of events compared to the expectation 
from the Standard Model in very short distance e+p deep-inelastic scattering pro- 
cesses has renewed the interest in the search for new physics which could manifest 
in electroweak-like interactions. New preliminary results from the HI and ZEUS 
experiments making use of all available e+p data are reviewed here, with an em- 
phasis on the search for new bosons possessing Yukawa couplings to lepton-quark 
pairs. The sensitivity of HERA to leptoquarks, and to squarks of R-parity vio- 
lating supersymmetry, is confronted to existing indirect constraints from rare and 
forbidden semi-leptonic decays, atomic parity violation and neutrinoless double- 
beta decay, as well as to direct constraints from LEP and Tevatron colliders. The 
HERA and Tevatron colliders are found to offer exciting prospects for new physics, 
accessing yet unexplored domains of the mass-coupling plane. Possible striking 
manifestation of explicit lepton flavour violation is also discussed. 

1 Introduction 

The interest in the search for a physics beyond the Standard Model (SM) 
which could interfere with electroweak-like interactions has been, considerably 
enhanced recently by the observation in the HI Id and ZEUS □ experiments 
at HERA of a deviation from SM expectations in neutral current (NC) deep- 
inelastic scattering (DIS) events at very high squared momentum-transfer 
and possibly, yet at a less significant level, in charged current (CC) processes. 
In particular, an apparent "clustering" of outstanding NC-like events at masses 
around 200 GeV haa. motivated considerable work on leptoquark constraints 
and phenomenology H and on squarks in_^p-SUSYu, while contact interaction^ 
were also discussed as a possible source of distortion of the high spectrad. 

The ep collider HERA offers unique possibilities to search for s-channel 
production of new scalar particles which couple to lepton-parton pairs, up to 
a kinematic limit of ^/s^ ~ 300 GeV. The mass range relevant for a possible 
discovert at HERA has now been severely constrained by the TeVatron exper- 
iments Uu in particular for scalars decaying with a large branching ratio Peg 
into electron+quark. These constraints can be partly avoided for the squarks 
in _^p-SUSY theories where Peg is naturally small given the competition with 
gauge decay modes. Exhaustive squark searches covering both^p decay modes 
and various possible gauge decay modes are thus strongly motivated. 
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The original HI and ZEUS results mentioned above were based on data 
samples collected from 1994 to 1996. The analyses have now been updated in- 
corporating also the 1997 data, for a total gain of a factor ~ 2.5 in integrated 
luminosity. New preliminary inclusive NC and CC HERA results using all 
available e'^p data will be reviewed in section |2| and |^ of this paper. Interpre- 
tation and constraints on contact interactions will be summarized in section ^ 
New results on leptoquark searches will be discussed in section || and searches 
for squarks of ^p-SUSY in section ^. Searches for direct lepton flavour violating 
processes will be discussed in the context of both leptoquark phenomenology 
and .^p-SUSY. The HERA results and prospects will be compared to those of 
indirect processes and to direct searches at other colliders. 



2 Nevif Physics and Very High Rates at HERA 

The investigations of very high DIS-like processes at HERA, even with 
low statistics, are strongly motivated by the potential reach for new physics 
beyond the SM such as the production of new scalar (Xs) resonances. An Xs 
produced in the s-channel could for example lead to individual event signatures 
indistinguishable from standard NC and CC DIS if it decays into e -I- q or 
1/ + q. The new signal would nevertheless be identified statistically as a peak 
in the invariant mass distribution, associated with a characteristic angular 
distribution of the decay products. As will be discussed in section ||, the Xs 
particles which decay uniformly in their CM frame would contribute most 
significantly at large or large y where = JVPy, M is the resonance 
mass and y is related the decay polar angle of the final state lepton. Such 
new bosons possessing Yukawa couplings to lepton-quark pairs can in addition 
contribute in the u-channel by converting a lepton into a quark (and vice versa). 
Interference with SM DIS processes could also originate from effective contact 
interaction caused by the exchange of new bosons of masses M ^ ^JSep■ 

It is with these considerationfi in mind that HI and ZEUS experiments have 
carried their original analysis Era of the 1994^96 data corresponding to inte- 
grated luminosities of ~ 14.2±0.3 pb^^ and CzEUS — 20.1±.5 pb^^. The 
detailed comparison of the high tail of NC and CC-like event rates with ex- 
pectation from SM DIS processes has revealed exciting features. These observa- 
tions have now been updated including 1997 data, for a total of Chi — 37 pb^^ 
and CzEUS — 46.6 pb""'^. I shall now briefly review the essential features of 
the original and updated rate measurements. The results obtained when "con- 
verting" the inclusive DIS measurements into differential cross-sections will be 
discussed in the next section. 

The selection of NC-like events at HERA is straightforward. It essentially 
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Figure 1: (a) QI distribution of the NC DIS candidate events for the 1994 -»97 HI 
data (•) and for standard NC DIS expectation (histogram); (b) ratio of the observed 
and expected number of events as a function of Q^; the lines above and below unity 
specify the ilcr levels determined using the combination of statistical and systematic 
errors of the NC DIS expectation; (c) and (d) : as (a) and (b) but for 1997 data alone. 



requires an isolated at high E^.e well balanced by the hadronic transverse 
flow and with no longitudinal losses "visible" along the e beam —z direction 
(.T,vis. E - 2^;^'='""). For any DIS-like event, as for the production of 

an Xs particle involved in a 2 ^ 2 body process, a mass M = ^SepX can be 
calculated. The Lorentz invariant a: represents (at lowest order) the momentum 
fraction of the proton carried by the incident struck quark. By kinematic 
constraints, M and y (hence Q^) can be reconstructed from two independent 
measurements such as the energy and angle of the final state lepton, or by 
combining two angles from the lepton and hadronic energy flowB. 

Fig. |l|a shows the distribution measured! by HI for the 1994— >97 data 
in comparison with the SM expectation for NC DIS. The ratio of the mea- 
sured over expected distributions is shown in Fig. ^3. The essential features 
observedtl in the original data remain visible. The NC DIS expectation is seen 
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Figure 2: Me and j/e distributions for HI NC DIS candidate events (•) at > 
15000 GeV^ for (a)(b) the 1994^96 data and (c) (d) the 1994^97 data. Superim- 
posed as open histograms are the standard NC DIS expectation. 
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to be in excellent agreement with the data for < 10000 GeV^. At larger 
Qg, deviations are observed, with a slight deficit in the range ~ 10000 — 
15000 GeV^ and a number of observed events in excess at Q^>15000 GeV'^. 
Similar plots are now shown in Fig. |l]c,d for HI 1997 data alone. It is seen 
that the new data "suggests" similar deviations but with marginal ^gnifi- 
cance. Considering the full set of e+p data, HI and ZEUS are left nx3 with 
shght excesses at highest Q^. For > 15000 GeV^, Hl+ZEUS observe 
7V^]+^^^'^ = 42 events while Ndis = 32 ± 8.5 events are expected. For 
Q2 > 20000 GeV^ A^o^]+^^^^ = 18 while Ndis ^ 9.5 ± 1. One expects an 
equal or larger upward fluctuation in 1% of random experiments. ZEUS most 
significant deviation comes from 2 outstanding events at > 35000 GeV^ 
from their original datasetH, where they now expect 0.29 ± 0.02. 

Fig. ||a and b show the HI and distributions at Ql > 15000 GeV^ 
for the original 1994— >96 dataset. The data was seen to exceed SM expectation 
most prominently around Me ~ 200 GeV for large ye- Out of the 12 HI events 
at Q2 > 15000 GeV^ 7 appeared to be "clustered" in the bin 200 GeV ± 
AM/2 with AM = 25 GeV where one would expect 0.95 ±0.2. This particular 
clustering was not specifically supported by the original ZEUS observationsllil. 
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Figure 3: Same distributions as in Fig. |2| taken from ref. □ but showing in addition 
to NC DIS, as an illustration, the expectation for a M = 200 GeV scalar leptoquark. 



I will come back on this question in section |5[ The same representation of 
HI Me and ye results but including 1997 data are shown in Fig. ^ and d. 
The excess in the mass bin at AI ~ 200 GeV is found to be less significant. 
Here the fulLdataset was re-analysed including a new in situe electron energy 
calibration E3I1j. This has led to slight migrations (within originally quoted 
systematic errors) of individual events in the M — y plane. In particular the 
Me values are measured on average to be 2.4% higher. Considering nevertheless 
the same central mass bin of 200 GeV ± AM/2 with AM = 25 where the most 
significant fluctuation was originally observed, HI now finds Nobs — 8 events 
while Ndis — 3.0 ± 0.5 are expected. Of these observed events, 5 originate 
from the 1994 96 data (for 38% of Chi) and 3 from the 1997 data (for 62% 
of Chi)- Hence, it is fair to say that the "clustering" around ^ 200 GeV 
is not confirmed by the 1997 data alone. 

Shown for illustration in Fig. ^ are Me and ?/e contributiona, expected for 
a scalar leptoquark of mass Mx = 200 GeV and given couplingEj. Such a new 
particle would only contribute significantly at highest Q^. Hence, in preserioe 
or not of such new physics, it is reassuring that HI and ZEUS have found uH 
an excellent agreement between data and SM expectation for the y or mass 
spectra when considering low minimal thresholds. Interestingly, at highest 
the shape of the original Me and ye distributions and the excess are seen to 
be very well described by a combination of NC DIS + leptoquark contributions. 

3 DiflFerential Cross-Sections at High at HERA 

The observed Qg rates can be turned into differential cross-sections by con- 
verting the measured number of events to bin averaged values, using Monte 
Carlo acceptance calculations and detector efficiencies. It should be empha- 
sized that such procedure implicitly assumes that the simulation of the SM 
model properly accounts for migrations (e.g. due to initial or final state gluon 
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Figure 4: Differential cross- 
section da/dQ^ for NC and 
CC DIS processes measured 
by the HI (squares) and 
ZEUS (points) with e+p data 
and compared to SM expec- 
tation (using here CTEQ4 for 
the parton momentum distri- 
butions). 
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radiation) from true to measured kinematic quantities. This assumption could 
be invalidated in ncesence of new physics such as resonant production of long- 
lived leptomesons The differential cross-sections da/dQ'^ extracted by the 
HIEI and ZEUS experiments are shown for NC and CC DIS in Fig. |[ 

Both experiments find for NC a remarkable agreement with SM expecta- 
tion for < 10000 GeV^ over 4 orders of magnitude. It should be recalled for 
completeness that in this range and at large x, the SM contribution for NC 
DIS proceeds dominantly through t-channel 7-exchange with a u valence quark 
of the proton. For CC DIS, due to the W mass term in the exchanged boson 
propagator, the cross-section is suppressed at low Q^, and falls less steeply 
than for NC in the intermediate range. In e+p, CC proceeds dominantly 
through t-channel exchange with a d valence quark of the proton. Thus, the 
CC cross-section remains a factor >4 below that of NC for Q^^IOOOO GeV^. 

Fig. ||a and b show for HI and ZEUS the ratio of the measured over ex- 
pected NC da/dQ'^ cross-sections. Here for the SM prediction, HI performsuLa 
NLO QCD fit while ZEUS uses the CTEQ4 parton momentum distributionsllj. 
To constrain the iugh x domain, the HI fit combines the F2 structure function 
data from NMclla^and BCDMSO on both proton and deuteron targets. It 
furthermore makes use of low Q2 hi data (Q^ < 120 GeV^) from 1994 and 
1995-96 0. The measurement is clearly seen here to be in excellent agree- 
ment with SM expectation in the range 1000 < < 10000 GeV^ As 
was already seen from "unbiased" event rates in Fig. ^, the HI measurement 
slightly undershoots the SM expectation at ~ 10000 GeV^ while an excess 
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Figure 5: Ratio of the mea- 
sured over the SM predicted 
NC da/dQ^ cross-section ob- 
tained by (a) HI and (b) 
ZEUS. Here for the SM, HI 
takes a NLO QCD fit extrap- 
olated from Q2 < 120 GeV^ 
data while ZEUS makes use 
of the CTEQ4D QCD evolved 
parton densities. The gray 
band shows (a) for HI the un- 
certainty on the absolute nor- 
malisation (i.e. measured C) 
and (b) for ZEUS the uncer- 
tainty on parton density func- 
tions. Inner (outer) bars cor- 
respond to statistical (total) 
errors. 
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is observed at the 2a level for Q2>i5nnn GeV^ For CC QI§7-|W and ZEUS 
had observed in their original analysis Era and still observe E3ii3^t2l a tendency 
for the data to lie above SM expectation at highest but the CC results 
receive large systematic error contributions from the experiments (hadronic 
energy scale) and "model" {d quark momentum density). For a measured 
Q\ > 15000 GeV^ HI observes 9 events for 5.1 ± 2.8 expected while ZEUS 
observes 8 events for 3.9l]^'g expected. 

The cPa/dxdQ^ cross-section extracted by HI is plotted as a function of 
for a wide range of fix x values in Fig. ^. The results are expressed in 
terms of the reduced cross-section aNC defined as o-nc = (a;Q^/27ra^) 1/(1 + 
(1 — j/)^) d^a/dxdQ^. They are compared to the MRST structure function fit 
extrapolated to high as well as to a NLO QCD fitB taking into account the 
HI high data. The high data is seen to further pull the NLO QCD fit 
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Figure 6: Measured HI /dxdQ^ reduced NC cross-section (•) compared to the 
MRST parametrization (upper curves) and to a new NLO QCD fit (lower curves) 
combining HI preliminary results with BCDMS and NMC. 



downward at highest relative to MRST. The upward "fluctuation" discussed 
above is_here consistently, visible at a; ~ 0.45 and Q^^ISOOO GeV^. From an 
analysiscj of the da/dx cross-section for > 10000 GeV^, and for the bulk 
of the cross-section which sits at the intermediate x range of x ~ 0.2 — 0.3, the 
suppression due to the negative 7 — interference in e'^p collisions (strikingly 
manifest from the inflexion of the theoretical pre dict i o n in Fig. ^ has now 
been solidly confirmed by both HERA experiments E£lE3. 



4 Searches for Contact Interactions 

Through the interference with SM gauge boson exchange, new bosons of mass 
M ^ JSep could affect the dajdQ^ cross-section measurements at HERA. 



Such a new interaction can be described as an effective 4-fermion "pointlike" 
{ee){qq) contact interaction (CI). CI was discussedEl as a possible "explanation" 
of an excess of very high events at HERA. The sensitivity to CI has now 
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Figure 7: (a) Ratio of the number of observed events per interval in ZEUS to 
the SM expectation. The full line shows the contact interaction (CI) best fit for the 
-|-AA model (see text); the shaded area is bounded by the ilcr limits, (b) Ratio of 
the measured da/dQ^ to the SM expectation for NC-Uke events in HI; the full line 
(dashed line) gives for a CI the 95% CL (resp. the CI best fit) for the S^j2- 



been investigatedtj using all available e^p data and considering the interaction 
Lagrangian dnt = Csm + ^c? where the chiral invariant C^f complement- 
ing the SM can be written in the form C^f = ^ {'nij{'^iliJ.^i){'i3l^'ij)] '^i*^ 
q = u,d and i,j = L, R and the coupling coefficients defined as rjij = ig^/A^'^. 
The -I- (— ) sign indicate constructive (destructive) interference. The "mass 
scale" A^ is conventionally defined when setting constraints as that scale rel- 
evant for a "strong" coupling strength of g'^ = A-k. 

The chirality structure of the CI model can_he chosen to avoid the severe 
constraints coming from atomic parity violationEil (AP V) . These are cancelled 
in particular if for the quarks g, rf^^ + 77^^ — r/^^^ -t- "rfj^j^ = 0, as realized for 
instance in the VV , AA and VA models considered by HI and ZEUS, with 
the mixture VV ^ LL + LR + RL + RR, VA = LL- LR + RL - RR, and 
AA = LL — LR — RL + RR. The SU(2) invariance is assumed which imposes 
vIl = vIl and r/]^^ = rf^^- 
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Table 1: Constraints 
at 95% CL on contact 
interaction scale A from 
HERA and other col- 
liders in models with 
VV, AA or VA structure 
(see text) for constructive 
(A"*") or destructive (A~) 
interference. 
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95% CL Lower Limits 
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Figure 0a shows CI constraints and best fit to the spectrum obtained 
by ZEUS for the AA(A+) scenario. Figure 0b shows as an iUustration the 
constraints and best fit obtained by HI in a leptoquark S'^j scenario. The 
•^1/2 (^^^ nomenclature in next section) possesses a —5/3 charge state coupling 
to e^u (e\u) and a —2/3 charge state coupling to e^d (ejd). Assuming 
SU(2) invariance, the CI couplings are t]\^ = 77^^ = -1/2(A/Afs)^ Both 
experiment rightly conclude E3 independently that no significant indication 
of a CI was found for these and various other models considered E3. It is 
nevertheless interesting to note that the best fits in the ^^(A+) scenario are 
found, consistently, to allow for a non- vanishing coupling. ZEUS finds ry/47r = 
1/A§ = 0.16l;^^ while HI finds 0.15^:^7 ; corresponding roughly to a value of 
Aq ~ 2.5 TeV. In practice, such a simple CI scenario is actually ruled out by 
TeVatron and LEP constraints as seen in Table 0. It is also interesting to note 
that, as seen in Fig. 0b an hypothetical scalar can be accommodated by HI data 
for Ms/ A = 387 GeV. For the VV and VA scenario, the HERA sensitivity 
is found to be comparable to that of other colliders. The VV{t^~) scenario 
would lead to a suppressicm of the NC dcr/dQ2 to a few 10^ GeV^ followed 
by a strong enhancement l3 qualitatively similar to the actual observation. 

5 Searches for Leptoquarks and Lepton Flavour Violation 

Generic leptoquarks (LQ) are colour triplet bosons which appear naturally in 
various unifying theories beyond the SM such as Grand Unified Theories and 
Superstring inspired i?g models, and in some Compositeness and Technicolour 
models. LQ searches have been carried either in the stcict context of the 
original Buchmiiller-Riickl-Wyler (BRW) effective modeled where the decay 
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Table 2: Scalar leptoquarks isospin families in the Buchmiiller-Riickl-Wyler model. 
These LQ will be in the following indexed with the chirality of the incoming electron 
which could allow their production, e.g. the Sq will be denoted by 5*0^. 



branching ratios are fixed by the model, or in the context of generic models 
allowing for arbitrary branching ratios. 

The BRW model considers all possible scalar (5/) and vector (V7 ) LQs of 
weak isospin /with dimensionless Yukawa couplings to lepton-quark pairs, 
where i and j indices denote lepton and quark generations respectively and L 
or R is the chirality of the lepton. The general effective Lagrangian obeys the 
symmetries of the SM and introduces 10 different LQ isospin multiplets, among 
which 5 are scalar families. These are listed in Table ^ in the so-called Aachen 
nomenclature and classification schemecj. The LQ search can be restricted to 
pure chiral couplings of the LQs given that deviations from lepton universality 
in helicity suppressed pseudoscalar meson decays have not been observed ES. 
This restriction to couphngs with either left- (A''') or right-handed (A^) leptons 
(i.e. A^ • A^ ^ 0), affects only two scalar LQs (S'o and 6*1/2). In all the results 
presented below, it is implicitly assumed as a simplifying assumptions that one 
of the LQ isospin doublet or triplet is produced dominantly and that the mass 
eigenstates within this multiplet are degenerate. 

5. 1 Search for first generation leptoquarks 

The search for first generation LQs at HERA involves the analysis of DIS- 
like events at very high Q^. The production cross-section glq depends on 
the quark momentum density in the proton and approximately scales with A^ . 
The scalar resonance which can then decay into e -f g or -I- g is expected to 
have a very narrow intrinsic width F = A^M/IGtt and the decays into a lepton 
and a quark jet lead to event signatures indistinguishable from SM DIS. A 
characteristic statistical signal of the direct production of LQs in the s-channel 
would be a peak in the reconstructed mass distribution. The LQ mass can be 
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reconstructed from the decay products (e.g. a charged lepton and a jet) or from 
kinematic constraints as M — ^Sf,pX. The isotropic decay of the LQs in their 
CM frame leads to a flat y spectrum. In this frame y = | (1 + cos 6**) with 9* 
being the decay polar angle of the final state lepton. This is markedly different 
from the da / dy ^ y~^ distribution expected at fixed x for the dominant t- 
channel 7-exchange in NC DIS events. Hence, the signal of first generation 
LQs would be most prominent at high y or equivalently at large = M^y. 

The HI LQ search for NC-like events requires an isolated e+ with a trans- 
verse energy of Et.e > 15 GeV. The reconstruction of the kinematics relies 
essentially (except for about ~ 4.4% of azimuthal range) on the e"*" energy and 
angle. The mass and y values of all events with > 2500 GeV^ are shown 
in Fig. |^(top left). The corresponding measured mass spectrum shape is found 
to be very well predicted by the SM as can be seen in Fig. ^(top right). Also 
shown there are the mass spectra after having applied a mass dependent y 
cut □ designed via Monte Carlo studies to optimize the signal significance for 
scalar LQ searches. This cut varies from ~ 0.6 at — 60 GeV to ~ 0.4 at 
Me ~ 200 GeV, and down to ye ~ 0.2 at Afe~250 GeV. HI observes 312 events 
satisfying this y cut in the mass range Me > 62.5 GeV, in excellent agreement 
with the SM expectation of 306 ± 23. In the mass range 200 GeV ± AM/2 
with AM — 25 GeV, Nobs — 8 events are found for Njjis ~ 3.0 ± 0.5. As was 
discussed in section ||, this slight excess mostly originates from the 1994-^96 
dataset. These 8 events have an average mass of Me ~ 206 GeV. When dis- 
cussing sensitivity and constraints on LQ searches, it should be kept in mind 
that according to LQ Monte Carlo simulations, the mass Me as deduced from 
the positron tends to systematically underestimate a true LQ mass by 2% 
for Mlq of 0(200 GeV). The most significant fluctuation is observed in HI 
for Me values which would correspond to Mlq ~ 210 GeV. 

The ZEUS LQ search for NC-like events requires a total Et > 60 GeV. 
A mass Mej is attributed to each candidate event and calculated from the 
positron and the jet at highest Et, not correcting for the finite jet mass. 
The Alej and cos 6** for all events are given in Fig. ||( bottom Left). Interest- 
ingly, the outstanding high events already discussed in ref. □ here tend to 
cluster around Mej ~ 215 GeV. The projected Mej spectrum is shown in 
Fig. ||( bottom right) after having applied conservatively a fiducial cut (shaded 
area in Fig. ^ bottom left)) which removes a difficult region for electron energy 
measurements at the interface between two calorimeters, thus cutting away 2 
very high events. In the remaining acceptance and for Mgj > 200 GeV and 
cos 6** > 0.25, ZEUS observes 7 events in slight excess of the SM expectation 
of 4.3 events. Of these 7 events, 5 originate from the 1994^96 dataset. For 
Mej > 200 GeV and in the full cos 9* range, 68 events are observed for an ex- 
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pectation of 43_]^2 from SM. This excess is found by ZEUSE2I to be mostly due 
to events having a cos 9* spectrum shape compatible with SM DIS expectation. 

While interesting upward fluctuations of the number of observed events 
have been found by both HI and ZEUS for masses Af >200 GeV in their 
1994— >96 dataset, the 1997 data disappointingly offered no confirmation of 
an excess in the NC DIS-like channel which would possess the characteristics 
of lepton-quark resonant production. Hence, both experiments proceed to set 
mass dependent constraints on the production cross-section ctlq of first gen- 
eration scalar LQs, treating the deviations observed as statistical fluctuations. 
Model "independent" upper limits on aLQ were derived by ZEUS in ref.113 and 
are shown in Fig. 0. 
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Figure 9: Cross-section exclu- 
sion limits at 95% CL from 
ZEUS for generic scalar and 
vector leptoquarks. 
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Further interpretatioD_pf such constraints can be discussed in the frame- 
work of the BRW model c3 where (3eq values are fixed by the model for the 
various LQ types, as given e.g. for scalars in Table ^j. The exclusion upper 
limits on aLQ translated into mass dependent limits on An as obtained by HI 
are shown in Fig. |l^ for scalar LQs having fermionic number F = 0, i.e. which 
can be produced via a fusion between the and a u or d valence quark. The 
limits are given here for masses up to 275 GeV above which one needs to move 
away from a resonance-like search since the mass peak of = scalars be- 
comes severely distorted. Very similar results are obtained by ZEUS who also 
derived constraints for vector LQs 113. The e+p collisions naturally offers the 
best sensitivity to _F = LQs and the limits obtained represent an irmrove- 
ment by a factor ~ 3 compared to previously published HERA results E3. For 
a coupling of the electromagnetic strength X^/An — aEM (i.e. A ~ 0.3), such 
LQ are excluded at 95% CL up to 275 GeV. The highest sensitivity (hence 
most severe constraints) is obtained for since both charge states can be 
produced via a e+u and e^d fusion. Only e^u (e+d) fusion is possible for the 
s-channel production of Sy^ {'^1/2)- '^^'^ F = 2 LQs will be best probed with 
the forthcoming e^p data taking starting in 1998 at HERA. 

The mass range of interest for a possible discovery at HERA of a LQ of 
the strict BRW model has now been severely reduced by the TeVatron pp 
experiments l3'Q, where first generation scalar LQs with M < 242 GfV and 
(3eq = 1 are excluded (95% CL) independently of the A (see ref. E3). For 
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Figure 10: Exclusion lim- 
its at 95% CL on the 
Yukawa coupling A as 
a function of the lep- 
toquark mass for F — 
scalar leptoquarks de- 
scribed by the BRW 
model. Domains above 
the curves are excluded. 
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f3eq = 0.5, the excluded domain reaches M ~ 200 GcV. The mass range 



excluded by the JDO experiment alone is shown in Fig. 10, At LEP e+e^ 
collider, searches l3 for single LQs using data taken at centre of mass energies 
of 161 and 172 GeV gives for A ~ 0.3 a best resulting limit of 142 GeV which 
is not yet competitive with the sensitivity achieved at HERA . 

Moving away from the BRW model, one can consider from a phenomeno- 
logical point of view scalar LQs undergoing both NC and CC DIS-like decays 
{Peg X Pi^q' 7^ 0) with and /3^qi treated as free parameters. By gauge in^ 
variance this can only be relevant for LQ coupling to e^u. HI has derived □ 
such constraints. For = 90% and Peg = 10%, LQ masses below 210 GeV 
are found to be excluded at 95% CL for a coupling An ~ 0.3. This extends far 
beyond the domain excluded by TeVatron experimentsQD which for such small 
values of f3eq only exclude scalar LQ masses below ~ 110 GeV. Given a value 
of All, upper limits on a^g can be translated for narrow intrinsic widths in 
terms of mass dependent limits on the branching f3eg without making specific 
assumptions on the nature of the other decay modes. Results are shown in 
Fig. [ll| a and b for LQs produced via e^d and e+u fusion respectively. Despite 
the small An values considered, the exclusion domains arCj^een in Fig. ^ to 
extend beyond the region covered by the Z50 experiment El at TeVatron for 
small Peg even in the less favourable case of an LQ coupling to e+d. HERA 
rules out masses below ~ 210 GeV if An = 0.1 (Fig. |ll|b) or below ~ 255 GeV 
for a interaction of electromagnetic strength (see for instance Fig. [l^a). 

5.2 Search for leptoquarks with mixed lepton flavour couplings 

We have seen in the previous section that the mass range of interest for a dis- 
covery of lepton-quark resonances at HERA is severely constrained by TeVa- 
tron results for LQs coupling solely to first generation fermions. Hence, HI 
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Figure 11: Mass dependent 
exclusion limits at 95% CL 
on the branching ratio P^q = 
BR{LQ —> eq) for scalar lep- 
toquarks produced by (a) e^d 
and (b) e^u fusion. Two ex- 
clusion regions (light dotted 
grey) corresponding to A = 
0.1 and A = 0.05 are repre- 
sented. The D% limit is also 
shown as hatched region. 
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has investigated El further the case of LQs possessing An x Asj ^ 0. Such a 
LQ coupling to both first and third generation leptons would lead to striking 
lepton flavour violating (LFV) processes e+g — > LQ r^q' . 

The analysis for the LQ t + q decays covers the hadronic decays of the 
T as well as the decay r+ fi^v^Dr but only the former is actually used in 
the exclusion limits derivation. For the latter, the analysis requires simply an 
isolated track with transverse momentum Pt > 10 GeV linked to the primary 
vertex and a visible transverse momentum flow imbalance as measured solely 
via calorimetry of PT^^iss > 25 GeV. Only four /i-f jet events are found to 
satisfy these criteria. AlLfour are amongst the "outstanding" high Pt lepton 
events discussed in ref. £3. None of these ii + X events has a /i candidate at 
9^ and i/)^ angles corresponding to the r angles predicted from kinematical 
constraints (the /i from the t decay should be strongly boosted in the r direc- 
tion). Hence neither LQ fj, + q (2"'' generation) nor LQ t + q; t ^ ^vv 
(3'"'' generation, /i channel) candidates are found. For the hadronic decays 
of the T, the analysis requires a "pencil-like" jet defined by a small invari- 
ant mass Mjet < 7 GeV and a low multiplicity of associated charged tracks 
(Ntracks < 2). The r jet furthermore must have either a large intrinsic tt*^ 
component or a stiff leading track, and be back-to-back in azimuth with the 
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Figure 12: (a) Mlq dependent exclusion limits at 95% CL on the Yukawa coupling 
All, for scalar leptoquarks produced by e^u fusion. Different hypothesis for the 
branching ratios into eq, rq axe considered. Domains above the curves are excluded, 
(b) Exclusion domains in the plane Asj {j — 1, 2) against Mlq for several fixed value 
of All (greyed areas). Indirect limits on A31 are represented by the dotted curve. 



rest of the hadronic flow. No candidate is found for these r hadronic decay 
modes while 0.8 ± 0.3 misidentified background event is expected. 

Mass dependent exclusion hmits are shown for An in Fig. |l2|a when fixing 
the relative Peg and (3rq branching fractions for a generic LQ couphng to e"*" + u 
pairs (such as the S^^^ l BRW model) and for three different sets of 

iPeq^Prq)- Here both channels LQ e"*" + jet and LQ r+ + jet are 
combined. This latter channel is essentially background free but the former 
benefits from a higher selection efficiency, such that both provide comparable 
sensitivity. Hence, when /3e<j + (irq — 1, the exclusion domain is very similar 
to the one obtained for P^q — 1- Assuming f3eq = 10% and Prq — 90%, masses 
below 275 GeV are excluded at 95% CL for An — V47ra as shown in Fig^ 12|a. 
An alternative representation of the LFV constraints is given in Fig. IIP" m 
the plane X^j against Mlq for different fixed values of An. Here again, a 
LQ formed via e'^u fusion is considered such that only couplings X^j with 
j = 1,2 are relevant. Furthermore the simplifying assumption Peq + Prq = 
1 is made. Fig. [T^ b shows the domains excluded at 95% CL by the r + 
jet final states analysis. Of course, for each hypothetical value of An, part 
of the mass and An x X^j domain is already implicitly excluded anyway by 
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the e + q analysis alone. Nevertheless, the r + jet analysis allows extend 
considerably the coverage of possible An x X^j values; e.g. for An — 0.03 
(1/100 of electromagnetic strength) the mass range above Mlq ~ 80 GeV 
is covered in complement. For An — Xsj ~ 0.03, these new HERA results 
extends the mass limits from previous HERA resultsEj by ~ 65 GeV. 

The best indirect constraintEl on A31, shown in Fig. [T^ b in the case An — 
0.3, is found to be typically one order of magnitude less stringent than HERA 
exclusion limit. This constraint comes from the upper limit on the branching 
ratio Pr^TT^e which could be affected by the process r — > d+LQ*; LQ* e + d. 
No low energy process constrains the coupling A32. Thus, a domain which was 
yet unexplored is now covered by HI in the plane A32 against Mlq. It should be 
noted however that other LQ species might suffer more severely from indirect 
constraints. This is the case in particular of LQs coupling to e+d pairs (such 
as the >§i/2,L in BRW model) for which the the couplings A31, Ago or A33 are 
constrained respectively by r ^ 7r°e, r K'^e and B —> reXcB. Such LQs 
have characteristics similar to the squark in SUSY models where R-parity 
is violated (see next section) by couplings X[ji (analogous to An) and Ag^^, 
(analogous to X^k)- Exclusion domains which have been derived inE3 in the 
plane Ag^^ against Mu show that even in such a case, HERA sensitivity extends 
beyond that of low energy phenomena for masses Ms < 250 GeV. 

The CDF experiment has performed a search for third generation LQ 
looking at rrbb final states EJ, and excludes a scalar LQ with masses belciw 
99 GeV if Prb = 1 • A complementary search has been carried out by Z?0 l3, 
where the analysis of i^vbh final states leads to a lower mass limit of 94 GeV for 
/3i,b = 1- These constraints are less stringent than the ~ 110 GeV mass limit 
obtained from TeVatron searches for first generation LQs assuming (3eq = 10%. 
Thus, a situation with small f3eq and large f3rq is especially favourable, leaving 
open e.g. for HERA an important discovery window at Mlq < ^Sep- 

6 Searches for Squarks of i?-parity Violating Supersymmetry 

Squarks are scalar SUSY partners of the quarks. The most general SUSY 
theory which preserves gauge invariance of the Standard Model (SM) allows 
for Yukawa couplings between one scalar squark {q) or slepton (I) and two 
known SM fermions. Such couplings induce violation of the i?-parity defined 
as Rp = (_i)3-B+L+2S^ where S denotes the spin, B the baryon number 
and L the lepton number of the particles. Of special interest for HERAEj are 
those Yukawa couplings A'^^j. (j, k are generation indices) violating the leptonic 
number and which couple to a sauark to a lepton-quark pair. The search for 
^p-SUSY at HERA was carried considering otherwise the field content of 
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Table 3: Squark production pro- 
cesses at HERA (e"*" beam) via a 
7?-parity violating A'^jj. coupling. 







pr 


oduction process 




111 


e+- 






e+H 


-ci- 


» ML 


112 


e+- 


f u 




e+H 


vs- 


UL 


113 


e+- 


f u 




e+- 


h5- 


ML 


121 




f c 


— > J* 


e+- 


Vd- 


> CL 


122 


e+- 


f c 


, S* 


e+- 


F s - 


> Cl 


123 


e+- 


f c 




e+- 


f 6- 


* Cl 


131 


e+- 


f ? 




e+ - 


f rf- 




132 


e+- 


f f 




e+ - 


f s - 




133 


e+ 


f i 




e+ - 


f fe- 





the Minimal Supersymmetric Standard Model (MSSM) and assuming that the 
neutralino is the lightest supersymmetric particle (LSP). 

Production and decay: With an e"*" in the initial state, each of the 
nine ^'ijk coupling allows for a resonant production of squarks through a e ~ q 
fusion process listed in Table ||. The squarks decay either via their A' coupling 
into SM fermions, or via their gauge couplings into a quark and a neutralino 
Xi [i = 1,4) or a chargino x'j U = 1>2). The mass eigenstates Xi and 
X^ are mixed states of gauginos and higgsinos and are in general unstable. 
This latter point holds in .^p-SUSY, in contrast to the strict MSSM, also for 
the LSP which decays via A'j^^j, into a quark, an antiquark and a lepton. In 
cases where both production and decay occur through a X[jf. coupling (e.g. 

Fig. ^3|a and c for ^ 0), the squarks behave as scalar LQs. The d^* can 
decay either into e+ + or + ■ Gauge invariance forbids the m^ vq 
decay and, hence, such squark type is left in the l^p decay mode with vPj^ — > 
+ d''. Hence, the final state signatures consist of a lepton and a jet and are, 
event- by- event, indistinguishable from the SM neutral and charged current 
DIS. In cases where the -u;^ (resp. d^*) undergoes a gauge decay into a Xa 
or a Xp (resp. x°), (e.g. Fig. and d) the final state will depend on the 
subsequent decay of the x- Neutralinos can undergo the^p decays Xa ~^ c^QQ' 
or Xa ~* the former (latter) being dominant if Xa dominated by its 
photino (zino) component. When x^ decays via l^p into a charged lepton, 
both the "right" and the "wrong" sign lepton (compared to incident beam) 
are equally probable leading to largely background free striking signatures for 
lepton number violation. On the contrary, the only l^p decays for charginos 
are xt ~^ vu^df' and Xa ~^ e^d'^d^ . Neutralinos Xa '^itih a > 1 as well as 
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Figure 13: Lowest order s-channel diagrams for first generation squark production 
at HERA followed by (a),(c) decays and (b),(d) gauge decays. In (b) and (d), 
the emerging neutralino or chargino might subsequently undergo 1^ decays of which 
examples are shown in the dashed boxes for (b) the a^nd (d) the Xi- 



charginos can also undergo gauge decays into a lighter x and two SM fermions, 
through a real or virtual gauge boson or sfermion. 

Event topologies: Overall, a e^+multijets final state configuration is a 
likely one for a significant part of the MSSM parameter space Y^^^fV^ considering 
squark gauge decays, as discussed in previous HERA analysisEjCj. Such event 
topologies have been searched by HI requiring an isolated found at large ye 
{ye > 0.4) while, nevertheless, being accompanied by 2 forward (relative to the 
incident proton direction) high Et jets. Such a configuration is very unlikely 
for standard NC DIS. HI finds 289 candidates in the e^+multijets channel, 
in good agreement with the mean SM background of 285.7 ± 28.0 expected 
from NC DIS and photoproduction (the latter contributes to less than 3%). 
The measured mass spectrum of the S3 selected events is compared to SM 
expectation in Fig. [ij. Good agreement is observed with only a slight excess 
of events observed at the highest masses. HI has further looked for "wrong" 
sign final state lepton in the process e^q' — > g — > e~q"q"q' . One such wrong 
sign candidate was observed in the data, while the SM prediction as estimated 
from Monte-Carlo was found to be 0.49 ± 0.2 (coming from NC DIS). Hence, 
apart from an excess of NC-like events at large Ql or large Me, no significant 
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Figure 14: Mass spectrum for e + mul- 
tijets final states for data (symbols) 
and NC DIS expectation (histogram). 
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deviation from the SM expectations has been found in a key topology for 
gauge decay channels. Assuming that the slight deviations observed are due 
to statistical fluctuations, the searches in IjLp and gauge decay channels can be 
combined to set constraints on ^p-SUSY models. 

Rejection limits: The rejection limits are derived as a function of the 
{(■^ mass assuming that only one of the ^'iji is non vanishing and combining 
all contributing channels. The masses of other sfermions are assumed to only 
influence weakly the branching ratios of the neutralinos and charginoso. Re- 
jection limits on as a function of squark mass are shown in Fig. |l^a for the 

v^j^ when combining the relevant event topologies, taking into account either 
NC-like e"''-|-jet only, or NC-like e^+jet combined with e^+multijets, or all 
three channels (including "wrong sign"e~-|-multijets). The MSSM parameters 
have been set here to ^ = -200 GeV, M2 = 70 GeV and tan/3 = 1.5. With 
this choice of parameters, the lightest neutralino is mainly dominated by 
its photino (7) component and M^o ~ 40 GeV, while the xt ^^id are nearly 
degenerate around 90 GeV. Combining the three contributing channels im- 
proves the sensitivity on ^'iji by up to a factor ~ 5 at lowest mass compared 
to the one obtained using only the NC-like channel. The relative contributions 
of the three channels in the case where Xi is 7-like are plotted against the 
squark mass in Fig. ^b for A' at the current sensivity limit. It is seen that 
for masses up to ~ 230 GeV, the channels e+-|-multijets and e~-|-multijets 
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Figure 15: (a) Exclusion up- 
per limits at 95% CL for A'lji 
as a function of the Mq, for 
a set of MSSM parameters 
leading to a 40 GeV Xi dom- 
inated by its 7 component. 
The limits are given for dif- 
ferent possible combinations 
of the contributing channels. 
Regions above the curves are 
excluded. (b) The rela- 
tive contributions of channels 
e++jet (SI), e+-|-multijets 
(S3) and e~-|-multijets (S4) 
versus Mq. 



have equal and dominant contributions. These channels play a decreasing role 
with increasing Mq as squark decays into Xi ^-^d X2 become kinematically al- 
lowed. The xi. and X2 become dominated respectively by their wino and zino 
components r4p^ and decay preferentially into vqq (a channel not covered in 
refc3 except partly through CC-like analysis). In the very high mass domain, 
a large Yukawa coupling is necessary to allow squark production, hence the 
relative contribution of e++jet is largely enhanced. Another set of values for 
(/.t, M2, tan/3) leadmg to a 40 GeV Xi dominated by its zino {Z) component 
was considered int3 to study the dependence of the rejection limits on the 
choice of MSSM parameters. In such a case the Xi of the Xi decay preferably 
in vqq (rather than in eqq'), leadingtojnultijets-l-P^'^j^g topologies not easily 
separable from the SM background Ej'ES and, hence, not expected to contribute 
very much to the sensitivity to new physics. Since the gauge decay width of 
the squark does not depend on the Yukawa coupling X'lji, the region of the 
plane (/3i, Mq) above the dotted line in Fig. |l5|b is excluded at 95% CL by HI 
combined analysis. In particular the branching ratio of a 200 GeV uj^ squark 
into e'^ + q is constrained to be smaller than ~ 1.5% (Fig. ^5|b) for the MSSM 
parameter choice presented here. It should be noted however that other spe- 
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Figure 16: Exclusion upper lim- 
its at 95% CL for the coupling 
X'lji as a function of squark 
mass, for various masses and 
mixtures of the Xii a-lso repre- 
sented are the most stringent in- 
direct limits on A'm and 
i = 2,3. 
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cific choices of (M2, /J,, tan/3) can allow for squarks at ~ 200 GeV to lead to 
NC-like topologies with /3e9~10%. 

Rejection limits obtained at HERA depending on the Xi niass and nature 
are compared in Fig. The sensitivity to for Mq < 200 GeV is better 

by a factor ~ 2 for a 7-like Xi than for a Xi dominated by its Z componepi, 
due to the highest part of total branching actually "seen" ia the HI analysisE3. 
One can infer from previous ^p-SUSY searches at HERAE3 that the two cases 
presented here are somewhat "extreme" and in that sense quite representative 
of the sensitivity at HERA for any other choice of MSSM parameters leading 
to a ~ 40 GeV Xi- The sensitivity to X[j^ for 7-like Xi increases with M^j 
given the corresponding increase of efficiency for the e-f-multijets channels. 
For Aiji — y/Anaem, squark masses up to 262 GeV are excluded at 95% CL by 
this analysis, and up to 175 GeV for coupling strengths ~ O.Olaem- For low 
masses, these limits represenLan improvement of a factor ~ 3 compared to HI 
previously published results E£l. 

Other indirect and direct constraints: The rejection limits obtained 
at HERA are compared to the best indirect limits in Fig. The most strin- 
gent indirect-constraint comes from the non-observation of neutrinoless double 
beta-,decayE3 but only concerns X'm coupling. The most severe indirect lim- 
its c2l on couphngs X'121 and A'^g^, which could allow for the production of 
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squarks c and i respectively, come from APVcil. It is seen that the sensitivity 
at HERA is better or comparable to the most stringent constraints on X'121 
and X'l^i- For large M^o values, HERA limits improves the sensitivity on some 
A' coupling by a factor up to ~ 4. 

LQ-like searches imply stringent constraints on uj^ squark masses only if 
somehow f3eq can be made large. But as explained above, this is unlikely in^p- 
SUSY for (say) Mg ~ 200 GeV. Hence, LQ-like constraints from the Tevatron 
are easily evaded. The problem is that a small Peg is so natural in^p-SUSYa 
that, at first glance, only a minute portion of the MSSM parameter space is 
left if one would like to "explain" a NC-like signal at M ~ 200 GeV in e+p 
collisions via the production of a 'U/,-like squark (e.g. or t^). Actually, 
elegant solutions can be found as discussed below in the case of the i^. At 
the TeVatron, SUSY searches have been mainly carried in the framework of 
minimal Supergravity (SUGRA) which imposes mass relations between the 
sparticles and i?-parity conservation. Recently, iD0E3 also considered squark 
pair production leading in 7?!p-SUGRA to like-sign dielectron events accompa- 
nied by jets, and has ruled out Mg < 252 GeV (95 % CL) when assuming 
five degenerate squark flavours. From a similar analysis by CDFc3 restricted 
to A'j2i 7^ o'^^ can infer that a cross-section five times smaller would lead 
to a Mg limit of ~ 150 GeV depending on the gluino and x° niasses. CDF 
also considered separately the pair production of a light stop ti assuming a 
decay into cxi and excluded Mj < 130 GeV. To translate this constraint in 
one relevant for X[^f. ^ 0, it should be noted that in this latter case, -^p-decays 
of the i would dominate over loop decays into cXi- Moreover, ^p-decays would 
themselves be negligible compared to i ^ bxt decays as soon as this becomes 
allowed, i.e. if M{ii) > M{xt) and if the ii eigenstate possesses a sizeable 
admixture of ti. The subsequent decays of the Xi would then lead to final 
states similar to those studied by CDF for ii cxi- Thus, 130 — 150 GeV 
appears to be reasonable rough estimate of the TeVatron sensitivity to a light 
t for A'j^gj, 7^ 0. In summary, TeVatron and HERA sensitivities are competitive 
in^p-SUSY models with five degenerate squarks, but models predicting a light 
t are better constrained at HERA provided that X[^j is not too small. 

Now let's come back to the case of the tL- Sizeable branching ratios for 
both I^p and gauge decay modes for a t produced via A'^g^ are very difficult to 
realize as argued ineil when restricting to i decays into e + d and b + xt ■ 
particular, the conflicting requirements due to APVEj constraints (implying a 
lower bound for /3eg) and to direct searches at the TeVatron (implying an upper 
bound on P^q) are not easily accommodated. A special case occurs if there 
exists a very heavy xt {M^+ > Mj:) and a light b (Afg < Alf) such that the t is 
left with the decay modes i — *■ e+d and i bW'^ . This interesting possibility 
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was first discussed in as a way to "explain" pSimultaneously an excess in the 
NC-like channel and the striking observationEj of LFV-like events with high 
Prmiss containing a high Pt muon and jet(s). For a hghtest 6 mass ~ 100 GeV, 
simultaneous sizeable branching ratios /Sg+d a-nd Piy^+ become possible 113 for 
the i thus extending the discovery potential at HERA for ~ 200 — 250 GeV. 

7 Conclusions 

The recent observation of possible deviations from Standard Model expectation 
in electroweak-like processes at HERA has considerably revived the interest in 
new theories requiring bosons with Yukawa couplings to lepton-quark pairs. 
Collider experiments and low energy precision experiments have reached re- 
markable (and often comparable) sensitivity to such particles, providing new 
avenues for the manifestation of exciting new physics. 
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